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; Marcello Cacciato 1 *, Edo van Uitert 2 , Henk Hoekstra 1 



O 

(N 



(N 

o 

u 

43 
6 



> 
in 

cn 
o 



1 Leiden Observatory, Leiden University, Niels Bohrweg 2, NL-2333 CA Leiden, The Netherlands 

2 Argelander-Institut fur Astronomie, Auf dem Hiigel 71, 53121 Bonn, Germany 



ABSTRACT 

The clustering of galaxies and the matter distribution around them can be described 
using the halo model complemented with a realistic description of the way galaxies 
populate dark matter haloes. This has been used successfully to describe statistical 
properties of samples of galaxies at z < 0.2. Without adjusting any model parameters, 
we compare the predicted weak lensing signal induced by Luminous Red Galaxies to 
measurements from SDSS DR7 on much larger scales (up to ~ 90 ft^Mpc) and at 
higher redshift (z ~ 0.4). We find excellent agreement, suggesting that the model 
captures the main properties of the galaxy-dark matter connection. To extend the 
comparison to lenses at even higher redshifts we complement the SDSS data with 
shape measurements from the deeper RCS2, resulting in precise lensing measurements 
for lenses up to z ~ 0.6. These measurements are also well described using the same 
model. Considering solely these weak lensing measurements, we robustly assess that, 
up to z ~ 0.6, the number of central galaxies as a function of halo mass is well described 
by a log-normal distribution with scatter (J\ a gL c = 0.146 ± 0.011, in agreement with 
previous independent studies at lower redshift. Our results demonstrate the value of 
complementing the information about the properties of the (lens) galaxies provided 
by SDSS with deeper, high-quality imaging data. 
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1 INTRODUCTION 

Since the advent of large and homogeneous galaxy surveys, 
it has become possible to constrain the relation between the 
observed properties of galaxies and their host dark matter 
haloes with ever increasing precision, albeit in a statistical 
sense. In particular, studies of the ob served abundances and 
clustering propertie s of galaxies (e. 



Griffiths et a! | Il99fj ; iHudson et all Il99gf ). the galaxy- 
galaxy lensing signal is now detected routinely as a 
function of the properties of the lens galaxies, thanks to 
multi-wavelength d a ta becoming readily available (e.g. 
i Fischer etall l2000l: iMcKav et alj l200ll; llloekstra et all 
bool 120051: ISheldon et al.l |2004|; IMandelbaum et all 1200^ 
Hevmans et al.ll2006l;|Parker et al.ll2007l ; iMandelbaum et al.l 



20081 : Ivan Uitert et all 120111 : iChoi et all \20lj ). The sys- 
tematics involved in these meas urements have also been 
studi e d in great detail (see e.g. IMandelbaum et"al 1 120051 . 
2006, 2008). The main application remains the study of 
the galaxy-dark matter connection, such as measurements 
of scaling rela tions between h al o ma ss and baryonic prop 
erties 
l200fj: 



2006: Yang et al. 2008: Moster et alj|20ld: Guo et al.l 


2010|; 


Behroozi et al.ll20ld: 


Moster et al.ll2012l: iGuzzo et al.l 


2000; 


Norberg et all l200ll. 


20021: Zehavi et al.l 120051: Wang et al. 



20071 ) have played a crucial role in establishing this relation 



with increasing detail. 

Complementing these methods, the weak gravita- 
tional lensing signal around galaxies of different observed 
properties (galaxy-galaxy lensing) has emerged as an- 
other powerful technique to co nstrain this re l ation . 
Since the first detections (e.g. iBrainerd et alj 1 19961 : 



E-mail: cacciato@strw.lcidcnuniv.nl 



(e.g. iHoekstra et all 120051: 



Ivan Uitert et all 
the halo 



Cacciato et al.l 120091: 



201 ll : IChoi et al l |2012|). 



Mandelbaum et all 

iLeauthaud et al.l l20ld : 

constraints 



prope rties (e.g. 



Mandelbaum et al 



Mandelbaum et al 



20061 : 



.Hoekstra et all 
Limousin et al.l ~ 



2004; 



2007; 



l200Sl ; Ivan Uitert et al.l |2012|) and 



meas u rements of bias param eters (e.g. Hoekstra et al.ll200ll . 
120021 ; ISheldon et all 120041 '). More recently galaxy-galaxy 
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lensing has also been used as a cosmological probe in 
combination with gal axy a bund a nce and/or clu s tering 
measurements dMore et al.l |2013|; I Cacciato et alj 120131 ; 



iRozo et alj|20ld : IZu et alj|2012l ; iMandelbaum et al.ll2012l l. 

There is growing scientific interest in probing the cos- 
mic evolution of structure formation in the Universe, which 
is now becoming possible thanks to new and forthcoming 
galaxy surveys. For instance, one can perform statistically 
repres entative analyses up to z ~ 1 i n the near future (e.g . 
KiDS Ide Jong et all I2012I. VIPERS iMarchetti et al. 1 12012| . 
Pan-STARRS iKaiser et"al]|2002l . DES^, HSC0), and possi- 
bly up to z ~ 2 in a decade, e.g. with missions such as 
LSST0, and EuclicQ (|Laureiis et al.ll201ll ). 

Alongside the progress in observational capabilities, 
theoretical modelling has also improved substantially. Nu- 
merical simulations have proven important to investigate the 
link between galaxy-g alaxy lensing and the g alaxy-dark mat- 
ter connection(e.g. Tasitsiomi et al. 2004 |_ Limousin et alj 



120051 ; iNataraian et al.ll2007l ; lHavashi fc White! 120081 ). Fur- 
thermore, the observed abundance, clustering and lens- 
ing signal have been successfully explained using a sta- 
tistical description of the dark matter distribution in 
the Universe as provided by the halo model (see e.g. 
ICoorav fc Shethl |2002| ; Ivan den Bosch et al.l 120131 ) coupled 
to a realistic model that describes the way galaxies of 
different observable properties populate host haloes (see 



e.g. |Y ang ct al .ll20Q. 



l2006l ; lYang et al.ll2008l ) 



le properties populate nost naioes (see 
iilCoorav fc Milosavlievic1l2005l ; ICooravl 



In this study, we examine the modeling of the galaxy- 
galaxy lensing signal up to z ~ 0.6. To this end, we 
first compare a model that describes the statistical prop- 
erties of galaxies at low redshift (van den Bosch et al. 
2013, Cacciato et al. 2013) to existing galaxy-galaxy lensing 
data me asured around Luminous Red Galaxies at higher 
redshift (Mandelba um et al.l 12012V To extend the redshift 
range even furthe r, and to obtain higher precision measure- 
ments, we follow Ivan Uitert et al.l (|201lh and complement 
the ninth data release (hereafter DR9) of the Sloan Digital 
Sky Survey (hereafter SDSS) with ~ 450 square degrees of 
high-quality imaging data from the second generatio n Red- 
sequence Cluster Survey (RCS2. lGilbank et alj|201lh . 

This paper is organized as follows. We describe the ana- 
lytical model in §2, its application to Luminous Red Galax- 
ies in §3. We then describe the surveys and the strategy 
to extract the new lensing measurements in Results are 
presented in [JS] Conclusions are drawn and discussed in SJB] 

Throughout this paper, we adopt the most basic 
('vanilla') ACDM cosmological model. Such ACDM cos- 
mologies are described by 5 parameters: the energy den- 
sities (in terms of the critical density) of baryons, fib, 
and cold dark matter, fidm; the spectral index, n, and 
normalization, ug, of the initial power spectrum; and the 
Hubble parameter, h 70 = H /(70 kms" 1 Mpc -1 ). The 
flat geometry implies that J1a = 1 — fi m = 1 — fib — 
fidm- Throughout the paper, following the results of Cac- 
ciato et al. (2013), we assume (fi m , fiA, o"s, hro, n, Q b h 2 ) = 



1 https://www.darkenergysurvey.org 

2 |http://ww.naoj.org/Projects/HSC/HSCProject.htmT1 

3 http://www.lsst.org 

4 http: / /www. euclid-ec.org 



(0.278, 0.722, 0.763, 1.056, 0.978, 0.0228). Radii and densities 
are in comoving unitfl. When physical units are used they 
are explicitly indicated with 'p— '. Furthermore, log is used 
to refer to the 10-based logarithm. 



2 MODELLING GALAXY-GALAXY LENSING 

In this section we briefly describe how model predictions 
for the galaxy-galaxy (hereafter g-g) lensing signal can be 
provided once one has a statistical description of dark mat- 
ter properties (i.e. their average density profile, their abun- 
dance, and their large scale bias) complemented with a sta- 
tistical description of the way galaxies of a given luminos- 
ity populate dark matter haloes of different masses (also 
known as halo occupati on statistics). The model is identical 
to the one presented in Ivan den Bosch et al.l J20131 ) and suc- 
cessfully applied to SDSS in Cacciato et alj (|2013l . hereafter 
C13). Readers familiar with this model may skip this section 
and continue from [J3] where we describe its application to 
Red Luminous Galaxies. 

Weak gravitational lensing is sensitive to the mass dis- 
tribution projected along the line-of-sight. Specifically, the 
quantity of interest is the excess surface density (ESD) pro- 
file, A£(i?), given by 



AE(i?, z lc ) = — I £(#, z lc ) R' dR' - £(#, z h 



(1) 



Here £(-R, z\ c ) is the projected surface mass density, which 
is related to the galaxy-dark matter cross correlation, 
£ g m(r, zi c ), according to 



E(R, Zi e ) = p n 



[1 + Cgm(r, zic)] dcj , 



(2) 



where the integral is along the line of sight with to the co- 
moving distance from the observer. The three-dimensional 
comoving distance r is related to uj through r = oj 1c + ui 2 — 
2wi e wcos#. Here, uji c is the comoving distance to the lens, 
and 6 is the angular separation between lens and source (see 
Fig.l in ICacciato et all [2009). Note that the galaxy-dark 
matter cross correlation is evaluated at the average redshift 
of the lens galaxies, z~i e . 

Observationally the ESD profile is inferred by measur- 
ing the average tangential distortion of background galaxies 
(sources) around foreground galaxies (lenses): 



(7t)(R) 



AE(i?) 



(3) 



where (...) indicates the azimuthal average inside an annulus 
at distance R from the centre of the lens and of width dR. 
In Eq.©, E cr it is a geometrical factor determined by the 
distances of (lens and source) galaxies: 



Ecrit 



A, 



47rG AcAc- SO (l + Zle) 2: 



(4) 



with Die, D so , and Di c _ ao the angular diameter distance to 
the lens, the source, and between the lens and the source, 
respectively, and the factor (1 + zi c )~ 2 accounts for our use 
of comoving units. 

5 We write the mean density of the Universe as p m = fimPcrit- 
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Under the assumption that each galaxy resides in a dark 
matter halo, AE(_R, z) can be computed using a statistical 
description of how galaxies are distributed over dark matte r 
haloes of different mass (see e.g. Ivan den Bosch et alj|2013h . 
Specifically, it is fairly straightforward to obtain the two- 
point correlation function, £gm(r, z), by Fourier transforming 
the galaxy-dark matter power-spectrum, P gln (k,z), i.e. 

£gm(r,z) = ^ J P gm (fc, z)^^k 2 dk, (5) 

with k the wavenumber. P gm (fc, z), can be expressed as a 
sum of a term that describes the small scales (one-halo, lh), 
and one that describes the large scales (two-halo, 2h), each 
of which can be further subdivided based upon the type of 
galaxies (central or satellite) that contribute to the power 
spectrum, i.e., 

P gm (fc) = P*(k) + P£(k) + P^(k) + P^ (k) . (6) 

As shown in lvan den Bosch et all (|2013l ). these terms can be 
written in compact form as 

P^(k, z)= tt x (fc, M, z) H y (k, M, z) n h (M, z) AM, (7) 



dMiHx(fc,Mi,z)nh(Mi,z) 



dM 2 H y (k, M a , z) n h (M 2 ,z) Q(k\Mi,M 2 , z) , (8) 



where 'x' and 'y' are either 'c' (for central), 's' (for satel- 
lite), or 'm' (for matter), Q(k\Mi, M2, z) describes the power 
spectrum of haloes of mass M\ and M2, and it contains the 
large scale bias of haloes as well as a treatment of halo 
exclusio n. Furthermore, ri h(M, z) is the halo mass func- 
tion of Tinker et al.l 120101 (see Ivan den Bosch et al.l 120131 ; 
ICacciato et al.ll2013l . for further detail). Here, we have de- 
fined 

M 

H m {k, M, z) = — u h (k\M, z) , 



H c (k,M,z)=H c (M,z) 



(Nc\M) 
n g {z) 



and 

K s (k,M,z)= u s (k\M,z) 
n g (z) 



(9) 
(10) 

(11) 



Here (N C \M) and (N S \M) are the average number of central 
and satellite galaxies in a halo of mass M = 47r(200p)P|oo/3, 
whereas n g (z) is the number density of galaxies at redshift 2. 
We compute these quantities using the following expressions: 



{N X \M) 



$ X (£|M) dL, 



(12) 



where $ X (L|A/) is the conditional luminosity function (see 
below and Appendix A), L_ and L+ refer to the lower and 
upper limit of a luminosity bin, respectively. Again, the sub- 
script 'x' stands for either 'c' (centrals) or 's' (satellites), and 



n g (z) = j {N e \M)n h (M,z)dM . (13) 
Furthermore, u s (k\M) is the Fourier transform of the 



normalized number density distribution of satellite galax- 
ies that reside in a halo of mass M, and u^{k\M) is the 
Fourier transform of the normalized density distribution of 
dark matter within a halo of mass M. In t his paper, sup- 
ported by the results of lCacciato et al.l i|2013l ) , we assume for 
both these profiles the functional form suggested in Navarro, 
Frenk & White (1997) . The conditional luminosity function 
($ X (Z/|M), hereafter CLF) describes the average number of 
galaxies with luminosities in the range L±dL/2 that r eside 
in a halo of mass M. Following ICacciato et all (|2013T ). we 
parametrize the CLF with nine parameters (see Appendix IA1 
for a thorough description). We note here that the CLF 
methodology describes the halo occupation statistics of both 
central and satellite galaxies and it is not limited to the 
choice of specific luminosity bins, rather it applies to galax- 
ies as a function of their luminosity. This will be of crucial 
importance when we will interpret the data presented in i|4] 



2.1 Additional lensing terms 

In the analytical model used by C13, which was summa- 
rized above, the lensing signal is modelled as the sum of 
four terms: two describing the small (sub-Mpc) scale signal 
mostly due to the dark matter density profile of haloes host- 
ing central and satellite galaxies; and the other two describ- 
ing the large (several Mpc) scale signal due to the clustering 
of dark matter haloes around central and satellite galaxies, 
respectively. This reads: 

AZ(R,z) = AEZ(R,z)+AT,lt(R,z) 
+ A^(R,z) + A^(R,z). 



(14) 



In the halo model the small scale signal (the 1-halo term) 
has two more contributors corresponding to: i) the baryonic 
mass of the galaxies themselves; and ii) the dark matter 
density profile of the sub-haloes which host satellite galaxies. 

The smallest scales probed by the data in this study 
are about 50 kpc, which are much larger than the typical 
extent of the baryonic content of a galaxy. Therefore, it is 
adequate to model the lensing signal due to the baryonic 
content of the galaxy as the lensing due to a point source 
of mass M g ~ Af star (see e.g. iLeauthaud et aUl2010T ). This 
reads 



AE lh ' s (7?,2) 



(M s tar(Zls))j 



(15) 



When accounting for the baryonic mass, this term adds to 
the other four indicated in eq. (|14p . Throughout the paper, 
we model the lensing signal as in Ej2j However, when de- 
scribing Figure O we comment on how model predictions 
are modified once the baryonic mass is taken into account 
in the simplified way described above. To that aim, we use 
the value of the average stellar masses, (M s t ar ) , for the galax- 
ies in the luminosity bins under investigation here (see §5). 
For completeness, we list these values in Table 1. 

The modeling of the dark matter density profile of 
the sub-haloes whi ch host satellite galaxies is conceptu- 
ally simple (see e.g . iMandelbaum et all 20051: Li et al.l 20091: 
i Giocoli et all l2O10l : Hi et al.ll2012l : iRodriguez-Puebla et all 
120131 ) However, a proper implementation of this term is ham- 
pered by the poor knowledge of the subhalo mass function 
(see e.g. Giocoli, Tormen & van den Bosch 2008) and of 
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the stripping mechanism (see e.g. iGao et al] |2004| ) which 
occurs once a dark matter halo enters a larger halo, i.e. 
when an initially central galaxy becomes a satellite. Many 
of the results about such subhalo properties are obtained 
from pure N-body simulations for which the limited mass 
resolution may still be a important limiting factor. Further- 
more, it is unclear how these results are affected by vari- 
ous baryonic processes in place during galaxy evolution (e.g. 
Ivan Daalen et aill201ll ). Given these uncertainties and since 
subhaloes only contribute a small fr action to the t otal lens- 
ing signal on small scales (see e.g. iLi et al.l |2009| ). in this 
paper, we refrain from modelling the lensing term due to 
the subhaloes which host satellite galaxies. We comment on 
the impact of this simplification when comparing model pre- 
dictions with actual measurements. 



3 SDSS LENSING SIGNAL AROUND LRGS 

The model summarized in §2 was used by C13 to fit the 
galaxy-galaxy lensing signal measurements performed via 
SDSS in the spatial range 0.05 < R < 2 Mpc and at red- 
shift z < 0.2. The same model can be used to make predic- 
tions about the scale and redshift dependence of the lensing 
signal. To test the robustness of the model, it is therefore 
interesting to examine how it performs, without any adjust- 
ments of the parameters (including the best-fit cosmology 
from C13), when compared to different data. 

We first consider the g-g lensi ng signal for a sampl e 
of Luminous Red Galaxi es (LRGs. lEisenstein et al.l l200lh . 
IMandelbaum et al.l (|2012h have measured the lensing signal 
around two LRG samples based on the SDSS DR7 catalogue. 
The selection of LRGs allows the study of the dark matter 
distribution via weak gravitational lensing at higher red- 
shift compared to the main sample. The effective redshifts 
of the two samples are z\ e ~ 0.26 and z\ c ~ 0.40. Both sam- 
ples have absolute magnitude limits —23.2 < M g < —21.1. 
Note that fc-corrections and evolution corre ctions to convert 
r-ban d magnitude to M g are taken from lEisenstein et alj 
(2001). Mor e details abou t the p rocedur e to select LRGs can 
be found in IKazin et all (|20ich and in IMan delbau m et al.l 
(|2012h . 

Figure [T] shows the data (filled circles with error bars). 
The model p redictions (solid lines) are obtained by using the 
same cuts as IMandelbaum et al.l (|2012h and the same model 
parameters found in C13. The model, although constrained 
using the main sample, describes the observed LRG signals 
very well. The value of the reduced \ 2 computed in the 
rang60 0.2 < R ~ 90/zf a Mpc is 1.0 for the main LRGs and 
0.8 for the high-z LRGs. The lensing signal around LRGs 
is reproduced over a large range of scales (0.2 < R < 90^ 
Mpc) and at high redshifts z ~ 0.26 and 0.4. 

It is worth emphasizing that the lensing signal on 
small and large scales carries different information. To first 
order, smaller scales probe the mass distribution within 
haloes, whereas larger scales probe the cosmological frame- 
work (mostly through a combinati on of the parameters 
Q m and erg). We recall here that ICacciato et alj (|2013l ') 



6 Note that the LRG lensing measurements used here are uncer- 
tain at scales smaller than about 0.14 Mpc (R. Mandelbaum 
private communication) 
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Figure 1. The e xcess surface density of LR Gs from the analysis 
of SDSS DR7 by IMandelbaum et alj d2012h . Solid lines refer to 
the model predictions for the lensing signal using the model pa- 
rameters retrieved in C13, without further adjustment. Note that 
the lensing measurements are uncertain at scales smaller than 
about 0.14 h^Q Mpc as indicated by the dotted vertical line. 



embedded their analysis in a fully Bayesian framework 
in which they also constrained the cosmological param- 
eters which define a 'vanilla' ACDM cosmology. They 
found that the parameters (fi m , ^a, f"8, , n, fl^h 2 ) = 
(0.278, 0.722, 0.763, 1.056, 0.978, 0.0228) best fit their model. 
Hence the agreement with the measurements is an impor- 
tant validation of the model determined by C13. It not only 
implies that the parameters that describe the halo occupa- 
tion distribution are also valid at higher redshifts, but also 
that the cosmological parameters are consistent. 

Before comparing the model to g-g lensing measure- 
ments based on a different data set in §4, we exploit the 
quality of the agreement between LRGs lensing data and 
model predictions to compute the average host halo mass of 
LRGs for both the main and the high-z sample. The estima- 
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Figure 2. Left: Redshift distributions of lens and source galaxies (with arbitrary normalization). The black solid histogram show the 
lens distribution of the DR9 sample used in this paper, whereas the blue dotted histogram refers to the lens distribution of DR7 used 
by van Uitert et al. (2011). The red dashed histogram indicates the approximate redshift distribution of the source galaxies. Right: 
Distribution of absolute magnitudes of the lens galaxies. The black solid histogram refers to the entire DR9 sample, whereas different 
lines refer to different subsamples defined via redshift cuts (see legend). 



tion of the average halo mass follows from 

<M 20 o) = —j—[{N c \M)n h (M,z lc )MdM, (16) 

n c (Z\e) J 

wher e rih(M, Z\ e ) is the halo mass function l| Tinker et alj 
2010) at the lens redshift, and (N C \M) is computed via 
eq. (|12[1 . We find that both the main and the high-z LRGs 
reside in haloes with (log[M2oo/(ftf 1 MQ)]) » 13.6, in gen- 
eral agreement with ind e pendent previous studies (see e.g. 
iMandelbaum et al.ll2006l : IZheng et al.1l2009l ). 



4 RCS2 GALAXY-GALAXY LENSING SIGNAL 

The weak lensing signal decreases as the lens galaxy ap- 
proaches the source. This limits the g-g lensing analysis of 
the SDSS main spectroscopic sample (DR7) to lenses with 
z ~ 0.2 when SDSS shape measurements are used. Reliable 
measurements at higher redshift are only possible by target- 
ing very luminous lens galaxies; see the case of LRGs in the 
previous section. A unique aspect of the SDSS is the wealth 
of spectroscopic information, which extends to higher red- 
shifts as is shown in Figure^ To improve the signal-to-noise 
ratio for higher redshift lenses, van Uitert et al. (2011) mea- 
sured the shapes of source galaxies using imaging data of 
a deeper survey that o verlap with the SDSS. Specifically, 
van Uitert et al] l|201ll) studied the lenses in the region of 
SDSS DR7 that overlaps with the second generatio n Red- 
sequence Cluster Survey (RCS2; lGilbank et al.|[201ll ). RCS2 
is a 900 square degree imaging survey in three bands (g' , 
r' and i') carried out with the Canada- France-Hawaii Tele- 



scope (CFHT) using the 1 square degree camera MegaCam. 
The RCS2 data are ~2 magnitudes deeper than the SDSS 
in the r'-band and the median seeing of 0.7" is roughly half 
that of the SDSS. Consequently the survey is well suited 
to improve the lensing constraints for lenses with z > 0.3, 
where the source distribution of the main sample in SDSS 
decreases significantly. In particular, it allows us to improve 
the lensing constraints for the most massive/luminous galax- 
ies, which are preferentially selected at higher redshifts. 

As shown in Figurc[2] a major change between the spec- 
troscopic samples of SDSS DR7 and DR9 is that high red- 
shift (z > 0.3) LRGs were targeted as par t of the Baryon 
Oscil lation Spectroscopic Survey (BOSS; lAnderson et al.l 
120121 ). As we show in this section, the lensing signal around 
these lenses can be determined with high precision using 
RCS2 shape measurements. 

The SDSS DR9 (lAhn et al.1 I2012T I overlaps with 471 
RCS2 pointings. This amounts to roughly 450 square de- 
grees (about 150 square degrees more overlap than between 
the RCS2 and the DR7 used in Van Uitert et al. 2011). 
The lens sample in the study presented here consists of all 
objects from the DR9 in the overlapping area that have a 
reliable spectroscopic redshift (according to the SciencePri- 
may flag) and th at are spectroscopically classified as galax- 
ies. In contrast to Ivan Uitert et aili|201ll ). we do not require 
the DR9 objects to have a match with an object from the 
RCS2 catalogues, but only that they reside within the field 
of view of a RCS2 pointing. This leads to a lens sample of 
~ 70, 000 objects, four tim es more than the lens sample used 
in Ivan Uitert et al.1 |201ll) . In the remainder of the paper, we 
shall refer to this sample as DR9. 
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Figure 3. Upper panel. The ESD signal measured in Van Uitert 
et al. (2011). Note the limited spatial scale (0.05 < R < 2 p-Mpc) 
and that the error budget is, on average, between 10 and 20 %. 
Lower panel. The ESD measurements presented in this study. 
Note that the probed spatial range now extends up to R ~ 10 
jj-Mpc and that the error budget is below 10 % for most of the 
probed scales. 



The redshift distribution of the lenses is shown in 
the left panel of Figure [2] For compariso n , we also show 
the lens sample used in Ivan Uitert et al] l|201ll) . labelled 
as DR7. The majority of lenses with redshifts z > 0.3 
are LRGs, whereas at lower redshifts, our lens sample 
consists of a mix of early-type and late-type galaxies. 



We do not consider these samples separately in this pa- 
per, as the halo model in use does not account for this 
split. The right-hand panel of Figure [2] shows the dis- 
tribution of absolute magnitudes for the whole sample, 
and for different redshift slices. The luminosities of the 
lenses are computed using the r-band Petrosian magni- 
tudes from the SDSS photometric c atalogues, corre c ted fo r 
extinction using the dust maps of ISchlegel et al.l (l998). 
K corrections were calculated to z = 0.1 using the 
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Figure 4. The ESD measured in this study for three subsamples 
of DR9: i) low (z < 0.2); ii) intermediate (0.2 < z < 0.5); and 
iii) high (z > 0.5) redshift. Note the high quality of the signal 
at both intermediate and high redshift, well beyond the regime 
probed by SDSS-alonc studies. 



KCOR RECT v4_2 code (jBlanton et al.ll2003l : lBlanton fc Roweisl 
Finally, a passive lu minosity evolution correc tion, 
E, was computed following iMandelbaum et al.l (|2012l ). In 
summary, the absolute magnitudes were computed as 
"^Mr = m + DM - K(z = 0.1) + E, where m is the appar- 
ent magnitude of a galaxy, DM is the distance modulus, 
K is the correction mentioned above, and E is the passiv e 
evolution correction taken from IMandelbaum et al.l |2012l ). 
i.e. E = 2(z - 0.1) [1 - (z- 0.1)]. We will comment on the 
impact of this assumption in §5.2.1 

The creation of th e shape measurement c atalogues for 
the RCS2 is detailed in Ivan Uitert et all l|201ll ) and we refer 
the reader to it for a detailed description. Since we lack red- 
shifts for the background galaxies, we select galaxies with 
22 < m r i < 24 that have a reliable shape estimat^ll as 
sources. The resulting average source density is 6.3/arcmin 2 . 
The approximate source redshift distribution for the sources 
(left-hand panel of Figure [2]) is obtained by applying iden- 
tical magnitude cuts to the photometric redshift catalogues 



7 We excluded sources with an ellipticity >1, and those whose 
photometry was deemed unreliable (e.g. due to image artefacts 
or neighbouring objects) as indicated by the flags of the source 
extraction program SExtractor (Betrin & Arnouts 1996). 
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of the COSMOS field (|llbert et al.ll2009f ). This procedure is 
detailed in Appen dix C. 

In contrast to Ivan Uitert et"al] (|201ll ). we do not limit 
our lensing measurements to individual pointings, but we in- 
clude the sources from neighbouring patches when measur- 
ing the azimuthally averaged tangential shear. This has the 
advantage that the lensing signal-to-noise at large radii im- 
proves, due to the larger number of sources at these separa- 
tions. Another advantage of including neighbouring patches 
is that it reduces the impact of systematic contributions to 
the lensing signal, as is explained in Appendix B. There we 
also present a detailed description of the error estimate. 

To compare the stat istical power between t he current 
work and the analysis in Ivan Uitert et all (120111 ). we show 
the lensing signal of the respective lens samples in Figure [3] 
We find that the signal-to-noise of the lensing measurements 
improves by about 50 per cent on average. Importantly, the 
lensing signal is robustly measured out to larger separa- 
tions (see Appendix B). We also split the lens sample in 
three redshift bins, and show the lensing signals for each 
bin in Figure [4] This illustrates that even at z > 0.5 (with 
zio ~ 0.59), we are able to obtain significant lensing measure- 
ments. Furthermore, the higher normalisation of the lensing 
signal measured at higher redshift is indicative of the fact 
that, not surprisingly, more massive lenses are selected at 
higher redshift. 



5 RESULTS 

The C13 model was constrained combining SDSS galaxy 
abundance and clustering measurements to g-g lensing 
data at low redshift, z < 0.2, and relatively small scales, 
R < 2Mpc. The comparison with the LRG sample in §3 pro- 
vides an important test of the model, but the lensing mea- 
surements are derived from the same pipeline as the data 
used by C13. It is therefore interesting to study how well 
the predictions compare to the results of the independent 
analysis that uses RCS2 shape measurements. Such a com- 
parison tests both the fidelity of the shape measurements 
and the model at even higher redshifts. To do so, we split 
the lens sample in eight luminosity bins and compare the g-g 
lensing signal to model predictions (see §5.1). Following this 
comparison we proceed to use the lensing measurements to 
investigate the possibility to constrain the galaxy-dark mat- 
ter connection at those higher redshift (see §5.2). 



5.1 Comparison of the RCS2 lensing signal to 
model predictions 

To study the model predictions as a function of luminos- 
ity we divide the DR9 lens sample in eight luminosity bins. 
The main properties of each bin are listed in Table 1. The 
black triangles with error bars in Figure [5] indicate the re- 
sulting ESD measurements based on the RCS2 lensing cat- 
alog as a function of projected lens-source separation. Over 
the large range in luminosity and scale (0.05 < R <, 10 h^ 
Mpc) we measure a significant lensing signal. Especially the 
measurements for the highest luminosity bins represent an 
improvement over what can be done with SDSS data alone. 
The magenta solid lines in Figure [5] correspond to the 
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Table 1. Properties of the excess surface density data 



Label 






lgM* 




•Nie 


(1) 


(2) 




(3) 


(4) 


(5) 


Li 


(-18.0, - 


17.0] 


9.21 


0.07 


1,418 


L 2 


(-19.0, - 


-18.0] 


9.72 


0.09 


3,650 


Ls 


(-20.0, - 


-19.0] 


10.23 


0.12 


8,918 




(-21.0, - 


-20.0] 


10.75 


0.19 


15,254 


L 5 


(-21.5, - 


21.0] 


11.09 


0.36 


14,013 


L e 


(-22.0, - 


21.5] 


11.32 


0.44 


13,555 


L- 


(-22.5, - 


-22.0] 


11.57 


0.51 


5,730 


Ls 


(-23.0, - 


-22.5] 


11.82 


0.59 


1517 



The galaxy samples used to measure the excess surface density 
profiles, AE(iJ). For each of these samples column (1) lists the 
magnitude label, column (2) lists the magnitude range, where 
M r = 0,1 M r , column (3) lists the log of the average stellar mass 
(lgM* = log [{Mstat/h^QMo}]), column (4) lists the mean red- 
shift, and column (5) lists the number of lens galaxies. 

model predictions based on the study by C13. Those predic- 
tions are in overall agreement with the RCS2 measurements 
of the ESD on all scales, including those well outside the 
range used to constrain the C13 model. The agreement be- 
tween model predictions and data supports the findings in 
C13 both in terms of the halo occupation statistics and the 
cosmological parameters. Furthermore, the mutual agree- 
ment of the model with t he data obtained with SDSS alone 
( Mandclb aum et al1l200rj ) and with those obtained here us- 
ing RCS2 implies overall consistency of the lensing signal 
from the two surveys. Interestingly, for the highest four lu- 
minosity bins, the C13 model predictions seem to system- 
atically over-estimate the g-g lensing signal at the la level. 
Although this aspect may be regarded as marginal given 
that the model predictions were not at all tuned to repro- 
duce these data, we will further investigate the relevance of 
this slight disagreement in §5.2. 

Model predictions based on C13 can be easily modified 
to account for the contribution to the lensing signal due to 
the galaxy baryonic mass (red dashed lines in Figure [5}. 
Using the simplest assumption for how the stellar content 
of the galaxy may contribute to the lensing signal (see fj3}, 
we note that fainter/less massive galaxies are virtually un- 
affected by such a correction, whereas brighter/more mas- 
sive galaxies exhibit a boost of the lensing signal on scales 
R <, 0.2/iyQ 1 Mpc, reaching up to a factor of about 1.5 on 
scales R ~ 0.05h^Q Mpc. 

We note here that in the current analysis the average 
stellar mass per luminosity bin is estimated by matching 
our lens sample to the MPA-JHU DR7 value added cata- 
logued which provides stellar mass estimates. Specifically, 
we use the matching objects to fit a linear relation between 
absolute magnitude and stellar mass, use this relation to as- 
sign a stellar mass to all our lenses, and finally determine 
the average for each luminosity bin (reported in Table 1). 
A technical caveat must be mentioned here: the MPA-JHU 
catalogue only contains galaxies from the DR7, and not the 
more recently observed ones from BOSS. Therefore, when we 
use the average stellar mass to compute the baryonic term 

8 http : //www . mpa-gar ching . mpg . de/SDSS/DR7/ 
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Figure 5. The excess surface density around lenses in DR9 that overlap with RCS2. The black triangles indicate the lensing signal 
measured using the RCS2 imaging data. Magenta solid and red dashed lines refer to the model predictions from C13, see ![2l and i|3] 
respectively. The shaded cyan region refers to the MCMC results where we consider Lqi^1i72 and o"i ogTjC as free parameters. 



in the halo model, we implicitly assume that the relation 
between luminosity and stellar mass is similar for the DR7 
galaxies as for those that were observed as part of BOSS. 
This assumption may not be accurate, but the use of the 
stellar mass in this paper serves only to roughly quantify 
on which scales and by what amount the g-g lensing signal 
might be affected by the baryonic mass of the galaxy. 

As a last cautionary note, we comment here on the fact 
that the model presented in this paper does not account for 
the mass distribution in the subhaloes which host satellite 
galaxies. As the quality of the lensing signal improves, it will 
become mandatory to add this extra term especially if one 



aims to retrieve the amount of stellar mass by fitting the 
small scale (R < 0.2hj Mpc) lensing signal. 



5.2 Constraining the galaxy-dark matter 
connection with weak lensing only 

More luminous galaxies reside on average in more massive 
haloes. Using Eq. (|16[) we find that this is indeed the case, 
and that the luminosity bins listed in Table 1 correspond to 
halo masses that range from (logM2oo/(/i7o 1 ^0) ~ 11-4, to 
(logAf20o/(^7o 1 ^o)) ~ 14-2. The signal-to-noise of the mea- 
surements presented in Figure [5] is highest for lens galax- 
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ies with —23 < 01 M r < — 21, which correspond to rela- 
tively massive haloes (12.5 < (logM 2 oo/(/iro M ©)> ~ 14-2). 
We explore here whether, thanks to the improved precision 
at the highest masses, we can constrain the model parame- 
ters which govern this regime using solely g-g lensing mea- 
surements. To this aim, we employ the same model used so 
far, but we now leave the parameters that govern the massive 
end of the galaxy luminosity-halo mass relation free to vary. 
In the parameterization used in this paper (see Appendix A 
for more detail), the relation between the luminosity of a 
central galaxy and its host halo mass is assumed to be: 



L C (M) 



L 



{M/Mip 



[l^iM/Kh)] 11 ' 12 
M N 72 



T ( M \ 



for M > Mi 



(17) 



Furthermore, the average number of central galaxies of a 
given luminosity is related to the halo mass via a log-normal 
distribution: 



$ C (L|M) dL 



log e 



2n a\ og L, 



-exp 



(logL - log L c ) 



(18) 



where 0~\ o& l c indicates the scatter in luminosity at fixed halo 
mass and log L c is, by definition, the expectation value for 
the logarithm of the luminosity of the central galaxy: 



logL c = / $c(L|M) logLdL. 



(19) 



Here, we consider Lo,M\, 72, and cri og L c as four free param- 
eters, while keeping 71 fixed to 3.18, the value retrieved in 
C13. The first free parameter has the units of a luminos- 
ity {K^qLq), the second has the units of a mass (A^Mq), 
whereas the remaining two are dimensionless. 

To determine the probability distribution of the model 
parameters discussed above we run a Markov Chain Monte 
Cark0 (hereafter MCMC) using the stand ard Metropolis- 
Hasting algorithm l|Metropolis et af1ll953[ ). In this chain, 
the parameters Lq , M\ , 72 and o-i og L c are f ree to vary and no 
prior information is used, whereas the remaining parame- 
ters are fixed at the same value as the one in the C13 model 
(see also Appendix A). As the satellite fraction is supposed 
to be very low for bright galaxies fe.g. [Mandelbaum et all 
l200fj : ICacciato et al.l2009l . 120131 : Ivan Uitert et al.l201ll ). and 
the faintest galaxies in this study have relatively large uncer- 
tainties, we do not expect significant biases from selecting a 
subsample of the model parameters that governs the galaxy- 
dark matter connection of central galaxies only. 

The 95% confidence levels of the g-g lensing models ex- 
plored with the MCMC are indicated by the cyan shaded re- 
gions in Figure[5] As expected, the subset of parameters that 
we have varied has almost no impact on the predictions for 
the lensing signal around the faintest galaxies. For brighter 
galaxies, the MCMC brings the model in better agreement 
with the observables than the initial C13 model predictions. 



9 The chain consists of four different chains which start from dif- 
ferent initial guesses in the parameter space. In total, we perform 
about three million model evaluations. With an average accep- 
tance rate of ~ 30%, the complete chain used in the analysis is a 
well converged chain of one million model evaluations. 




Figure 6. Posterior distribution of the parameter <7i og £ c . Blue 
shaded histogram refers to the analysis in §5.1 and §5.2, whereas 
the grey shaded histogram refers to the result of the test per- 
formed in §5.2.1 to assess the sensitivity of the analysis to passive 
evolution correction. The value of the scatter cri og £ c is robust 
against the uncertainties deriving from passive evolution correc- 
tion. 



The agreement between model predictions and data has im- 
proved by assigning smaller halo masses to galaxies of the 
same luminosity. From our analytical model (see especially 
eq. 17 and 18), one can see that lower halo masses at the 
same luminosity can be obtained by altering the L C (M) rela- 
tion at the massive end or by increasing the scatter, <7i Dg i c . 
As outcome of the MCMC we find that the L C (M) relation 
has substantially changed from the one retrieved in C13. 
However, as discussed in the following subsection, the infer- 
ence of the parameters which govern the L C (M) relation is 
very sensitive to the assumed correction for luminosity evo- 
lution, which is uncertain. Interestingly, the inference of the 
parameter o-i og L c is more robust against those uncertainties. 
Therefore, we report here only the corresponding result. The 
blue shaded histogram in Fig[6] shows the posterior distri- 
bution of the scatter in the number of galaxies of a given lu- 
minosity at any halo mass, cri og L c (see Appendix A for more 
details on this parameter). We find that a\ og i JC — 0.146T0.011 
(median ± one standard deviation), in excellent agreement 
with independent studies based on abundance, clustering, 
and/or satellite kinematics at lower redshift. Specifically, us- 
ing a large SDSS galaxy group catalogue, Yang, Mo & van 
den Bosch (2008) obtained ai ogLc = 0.13 ± 0.03 (black star) 
and they did not find evidence for a halo mass dependence. 
Cooray (2006) explicitly assumed no mass dependence in 
c"io g L c when studying the luminosity function and clustering 
properties of S DSS galaxie s , and found o-i og L c = 0.17ito'oi 
(grey triangle). iMore et al.l (|2009l ) studied the properties of 
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Figure 7. Upper panel. The reference absolute magnitude, 
01 M r = z M r + E, as a function of redshift computed via the pas- 
sive evolution correction, E, suggested by Blanton et al. (2003, 
blue solid line labelled B03) and Mandelbaum et al. (2012, dashed 
red line labelled M12) for the case of a galaxy with z M r = —20. 
Note that a galaxy with z M r = —20 at z > 0.1 will be rescaled 
to a fainter reference magnitude, 0,1 M r , at the reference redshift 
z = 0.1. Lower panel. The difference between the B03 and M12 
passive evolution corrections, A 0A M r = 0A Mf 03 - 0A M^ 112 . 



(2003) expression for the passive evolutio n would have only a 
minor impact on the lensing analysis by iMandelbaum et al.l 
as LRGs are selected only up to z ~ 0.4 and they are 
not split further into luminosity bins. However, in our anal- 
ysis lens galaxies are selected in narrow luminosity bins and 
up to 2 ~ 0.6. A different selection of lens galaxies directly 
translates into a different lensing signal. Spec ifically, we find 
that u sing the passive evolution correction of iBlanton et al.l 
(2003) leaves the lensing signal of faintest galaxies virtually 
unaltered, while leading to a higher normalisation (at about 
la level) of the lensing signal of the four brightest bins. 
Repeating the MCMC analysis as in §5.2 but on this new 
selection of lenses, we recover similar values of the scatter, 
<7 log l c , as indicated by the grey shaded histogram in Fig(6] 
However, we retrieve values for Mi, La, and 72 that differ 
by about 3 a from our initial analysis 10 !. The changes in 
the retrieved model parameters are to be attributed to the 
expected degeneracies due to the assumed parametrization 
(see eq. [17]). We conclude that the retrieved values of the 
model parameters Mi, Lo, and 72 are sensitive to the details 
of luminosity evolution correction. Therefore, in this paper 
we refrain from drawing any conclusion about the cosmic 
evolution of the L C (M) relation, deferring it to the analysis 
of a homogeneous sample of early-type galaxies (van Uitert 
et al., in preparation) for which passive luminosity correction 
can be modelled with higher confidence than in the current 
study where we consider a mix of early- and late-type lens 
galaxies. 



satellite galaxy kinematics around massive/luminous central 
galaxies and found o"i og L c = 0.16 ±0.04 (magenta circle). Fi- 
nally, C13 combining abundance, clustering and lensing of 
galaxies in SDSS found o"i og L c = 0.157 ±0.007 (red square). 

5.2.1 Sensitivity to passive evolution correction 

Ideally, one would like to compare the results on the L C (M) 
relation obtained here with those obtained at lower redshift 
to infer an evolutionary scenario. However, the physical in- 
terpretation of the results of the MCMC is hampered by 
the fact that the sample of galaxies used in this analysis is 
not uniform. In fact, the mix of early- and late- type galaxies 
changes with redshift due to the luminosity-based selection 
of the lenses. This leads to an uncertain correction for lumi- 
nosity evolution which enters in the definition of the refer- 
ence absolute magnitude of a galaxy, 0,1 M r (see §4). In what 
follows, we will show how a small variation in the correction 
of luminosity evolution impacts the lensing analysis and the 
corresponding model parameters. This example highlights 
the sensitivity of our model to the intrinsic luminosity evo- 
lution of a galaxy. 

In the analysis presented in §5.1 and 5.2, we have 
used the passive e volut ion correcti on, E, sugges t ed in 
IMandelbaum et al.l (|2012l ). However, IBlanton et ail (|2003l ) 
suggested a simpler functional form that has been widely 
adopted in the literature: E — 1.6(2 — 0.1). The difference 
between these two functions is highest at higher redshift, 
reaching about 0.3 magnitude at the highest redshift of in- 
terest here, z ~ 0.6 (see Fig. [7J. Using the IBlanton et all 



6 CONCLUSIONS 

We investigated how measurements of the galaxy-galaxy 
lensing signal around lenses at increasingly higher redshifts 
and at larger projected distances can be used to study the 
galaxy-dark matter connection. We showed that the an- 
alytical model presented in van den Bosch et al. (2013) 
and constrained by (SDSS) abundance, clustering, and lens- 
ing data at z < 0.2 (see C13), reproduces, without fur- 
ther adjustments, the galaxy-gala xy lensing signal measured 
around Luminous Red Galaxies (Mandelb aum et al.ll2012l ) 
at z ~ 0.26 and 2 ~ 0.4 and throughout the probed spatial 
range, 0.02 < R < 90 h^ Mpc (see Figure [J). This agree- 
ment is an important validation of the model determined by 
C13. It not only implies that the parameters that describe 
the halo occupation distribution are also valid at higher red- 
shifts, but also implies consistency with the cosmological 

parameters fo und by C13. 

Following lvan Uitert et all l|201ll ) , we measure the lens- 
ing signal around lenses from the Sloan Digital Sky Sur- 
vey (Data Release 9) using shape measurements from the 
450 square degrees that overlap with the second generatio n 
Red-sequence Cluster Survey (RCS2. Icilbank et alJlioiTI ). 
The higher source density and redshift results in a signifi- 
cant improvement, compared to SDSS data alone, for lenses 
with 2 > 0.3. We split the lenses into eight luminosity bins 
and measure robust tangential shear signals as a function of 
the transverse separation, R, in the range 0.05 < R < Wh^ 
Mpc (see Figure [5]). 



We have checked that also this MCMC has converged. 
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Compared to the earlier study by Ivan Uitert et al.1 

(2011) which used the overlap with DR7, the use of the 
overlap with DR9 increases the number of lenses, resulting 
in an improvement of about 50% in the precision of the lens- 
ing shear over the entire range probed here (0.05 <, R <, 10 
Mpc). In addition we now include the sources from neigh- 
bouring RCS2 pointings (previously the analysis was done 
on a pointing-by-pointing basis). This increases the lens- 
ing signal-to-noise at large projected lens-source separations 
(sec Figure [3}, and reduces systematic contributions to the 
lensing signal (see Figure iBl) . Finally, compared to the DR7 
catalogue, the redshift distribution of lens galaxies in DR9 
has a large number of galaxies at z > 0.4 (see Figure [2}, 
enabling us to probe the matter distribution at those high 
redshifts (see Figure 3}. 

We split the lens galaxies in 8 luminosity bins, rang- 
ing from -18 < 01 A/ r < -17 to -23 < 01 M r < -22.5. 
Brighter galaxies are distributed at increasingly higher red- 
shift such that the data span a wide range in redshift from 
z = 0.07 to z = 0.59. Moreover, since brighter galaxies live 
on average in more massive haloes, the range in luminosity 
probed here spans a correspondingly wide range in host halo 
mass. As a result, the measurements presented here simul- 
taneously probe the matter distribution in different regimes 
from small groups to massive clusters, and from low to high 
redshift (see Figure [SJ. 

Without any adjustment, the C13 model also describes 
the lensing signal obtained with RCS2 data very well. This 
corroborates the results based on the SDSS analysis of 
LRGs, but also implies consistency of the measurement 
of the lensing signal. We note that on the smallest scales 
probed here (0.05 < R < 0. 2h?J Mpc), we find better agree- 
ment if we include a contribution from the stellar mass of the 
galaxies: a simple point-mass model for the stellar compo- 
nent of the galaxies is sufficient to boost model predictions 
at those small scales leading to a better agreement with the 
data. 

Finally, exploiting the high signal-to-noise ratio of the 
lensing signal around bright galaxies, we attempt to con- 
strain aspects of the galaxy-dark matter connection across 
cosmic time. While the inference of an evolutionary scenario 
for the galaxy luminosity-halo mass relation is hampered by 
current uncertainties in the evolution of galaxy luminosity, 
we robustly assess that, up to z ~ 0.6, the number of central 
galaxies as a function of halo mass is well described by a log- 
normal distribution with scatter, 0i O gi, c = 0.146 ± 0.011, in 
agreement with previous independent studies at lower red- 
shift. 

Our results demonstrate the value of complementing the 
excellent information about the properties of the lenses pro- 
vided by the SDDS with deeper, high-quality imaging data. 
This allows us to probe the link between galaxies and matter 
around them in increasing level of detail and at increasingly 
higher redshift. In this paper we tested the model of C13 
and found that it overall performs very well. In future pub- 
lications we will use our data to examine the evolution of 
early-type galaxies only, and we will carry out a comprehen- 
sive study of the possible evolution with cosmic time of the 
galaxy luminosity-halo mass relation for early-type galaxies. 
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APPENDIX A: THE CONDITIONAL 
LUMINOSITY FUNCTION 

Throughout the paper, the average number of galaxies with 
luminosities in the range L ± dL/2 that reside in a halo of 
mass M is described by the conditional lu minosity function, 
$(L|M), introduced by I Yang et ail (120031 ): 



<A X |M> = 



$ x (L|M)dL. 



(Al) 



Li 



Follow ing ICoorav & Milosavlievid ([2005) and ICooravl 
(2006), we split the conditional luminosity function (here- 
after CLF) in two components, 



$(L|M) = $ C (L|M) + $ S (L|M) . 



(A2) 



where $ C (L|M) describes the contribution due to central 
galaxies (defined as those galaxies that reside at the cen- 
ter of their host halo), while $ S (L|M) characterizes satellite 
galaxies (those that orbit around a central). 

Our parameterization of the CLF mode l is motivated 
by the results obtained by lYang et al.l feOQSl ) from a large 
galaxy group catalogue (|Yang et al.l 120071) extract ed from 
the SDSS Data Release 4. and by iTal et all (|2012l ) from a 
study of the luminosity function of satellite galaxies of lumi- 
nous red galaxies. In particular, the CLF of central galaxies 
is modeled as a log-normal function: 



$ c (i|M) dL : 



log 



2"7T ai og L c 



-exp 



(logL - log L c ) 



dL 



(A3) 



and the satellite term as a modified Schechter function: 



$ 3 (L\M)dL = <t>; (i^)' 



exp 
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which decreases faster than a Schechter function at the 
bright end. Note that L c , a c , (j>g, a s and LI are in principle 
all functions o f the halo mass M. 

Follo wing Cacciato et al.l (|2009[). and mo t ivated by th e 
results of lYang et al.l (|2008h and iMore et ail l|2009l . l201ll ). 
we assume that o"io g L c , which expresses the scatter in log L 
of central galaxies at fixed halo mass, is a constant (i.e. is 
independent of halo mass and redshift). Note though that 
this does not imply that the scatter in halo mass at a fixed 
lumin osity, Q"i og M, is constant: as discussed in lCacciato et all 
( 2009) and IMore et~aH (2009), ciogM increases because the 
slope of the L C (M) relation becomes shallower with increas- 
ing M. In addition, for L c , we adopt the following parame- 
terization: 



L C (M) = L 



(M/Mi 



[1 + (M/Mi)] n 



(A5) 



Hence, L c oc M 71 for M < A/i and 
M 3> M\. Here Mi is a characteristic 
Lo = 2 71 ~ 72 L c (Mi) is a normalization, 
galaxies we adopt 

L* S {M) = 0.562L C (M) , 

Q S (M) = a s 



L c oc Af 72 for 
mass scale, and 
For the satellite 



(A6) 
(A7) 



(i.e., the faint-end slope of <!> s (L|Af ) is independent of mass 
and redshift), and 

log# s *(A/)] = 6o + 6i(logMi 2 ) + 6 2 (log Mwf , (A8) 

with M12 = M/(W 12 h^ M Q ). Note that neither of these 
functional forms has a physical motivation; they merely 
were found to ad equately describe the results obtained by 
I Yang et al.l l|2008l ) from the SDSS galaxy group catalogue. 

To summarize, our parameterization of the CLF thus 
has a total of nine f ree p arameters. Based on the re- 
sults of lCacciato et all (|2013l ). unless otherwise specified, we 
adopt the values (log[Mi/(/if 1 M Q )], \og\L /{h^ 2 L®)], ji, 
72, criog Lc , a s , b , 6i, 6 2 )=(11.39, 10.25, 3.18, 0.245, 0.157, 
-1.18, -1.17, 1.53,-0.217). 



APPENDIX B: SHEAR SYSTEMATICS AND 
COVARIANCE MATRIX 

For this paper the inaccuracies in the correction for PSF 
anisotropy is a main source of bias. On small scales we aver- 
age over many lens-source pairs which have a random orien- 
tation with respect to the direction of PSF anisotropy, and as 
a result the lensing signal is robust. At large radii the lensing 
signal is small and residual systematics may become more 
dominant because the angles between the PSF anisotropy 
and the lens-source pairs may no longer be isotropic because 
of masks or the survey geometry. In this section we there- 
fore examine the reliability of the lensing signal on scales 
> l/i^p-Mpc. 

To account for residual systematics in our shape mea- 
surement catalogues, and for the impact of image masks on 
the lensing signal, we compute the lensing signal around a 
large number of random points and subtract that from the 
galaxy-galaxy lensing signal. In the absence of systematics 
or an isotropic orientation of lens-source pairs, this signal 
should vanish. The red line in Figure iBll compares this sig- 
nal to the observed lensing signal (triangles with error bars) 



for the various luminosity bins. We find that the correction 
is generally very small: it is negligible for the LI to L5 bins, 
and for the other bins it is significantly smaller than the lens- 
ing signal over the range 0.05 < r < 10 h^ l p— Mpc which 
we use in our analysis. Therefore, any small error in the cal- 
culation of the random shear signal will have a minor effect 
at most on our results. 

The random shear signal is small because we in- 
clude neighbouring pointings in our lensing analysis. Con- 
sequently, at large projected separations, the lensing signal 
is averaged over many more lens-source orientations, which 
averages out any residual systematics in our shape measure- 
ment catalogues at those scales. That it is important to in- 
clude the neighbouring pointings, as opposed to analysing 
each pointing separately, is demonstrated by the dashed or- 
ange curves: these show the random signal for the case where 
we perform the lensing measurements on single exposures 
only. In this case, the random signal strongly increases with 
increasing lens-source separation, and its amplitude becomes 
comparable or even larger than the size of the lensing signal 
itself. Even a relatively small error in the computation of 
the random lensing signal could seriously affect the results, 
which is clearly undesirable. Therefore, it is important to 
conduct the analysis on patches rather than individual ex- 
posures. 

Another commonly used test is the measurement of the 
cross-shear, which is the projection of the source elliptic- 
ities to a unit vector that is rotated by 45 degrees from 
the vector that connects that lens and the source. Galaxy- 
galaxy lensing only produces tangential shear and not cross 
shear, as the average gravitational potential of a large num- 
ber of lenses with random orientations is circularly sym- 
metric. Therefore, if we measure a cross shear that is not 
consistent with zero this indicates the presence of systemat- 
ics in our shape measurement catalogues. For all luminosity 
bins we find that the cross-shear is consistent with zero on 
all scales used in our analysis. 

In order to quantify the level of correlation between the 
lensing signals of different radial bins, we compute the corre- 
lation matrix from the data using the "delete one jackknife" 
method (Shao 1986). We treat each pointing of the total N 
pointings as a sub-volume that is subsequently left out to 
create a new realization of the data. The covariance matrix 
is then determined as 



N - 



N 



1 

-^(7^-70(7* 



7i) 



(Bl) 



with ji the shear of the i-th radial bin, 7/ the shear at that 
location from one of the jackknife realizations, and 7i the 
mean shear of that bin determined by averaging over all the 
jackknife realizations. The correlation matrix follows from 

CoTTij = C-ij I y/ Cii Cjj • 

We find that the correlation matrix is practically di- 
agonal for almost all of our luminosity bins. Only for L2, 
L3 and L4 we find some low-level off-diagonal terms only 
around scales of ~1 h^p— Mpc. Since our analysis is mostly 
sensitive to the highest luminosity bins, we assume that the 
correlation matrices are diagonal when we fit the models to 
the data. 

Note that the correlation matrix that results from the 
jackknife method depends on the size of the sub-volume 
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Figure Bl. Assessment of systematics in the excess surface density measurements. Each panel refer to a luminosity bin as indicated 
by the labels. Black triangles with error bars refer to the excess surface density signal. The orange dashed lines denote the random shear 
signal obtained if one were to use only one single exposure, whereas the red solid line refers to the random shear signal obtained including 
neighbouring patches as in our analysis. 



that is subse q uently left out. This is demonstrated in 
iNorberg et al.l l|2008l ), who compared several ways to de- 
termine the variance and covariance of 2-point clustering 
measurements. For our purposes, the covariance matrix we 
determine is expected to be sufficiently accurate. However, 
for using measurements like these to constrain cosmologi- 
cal parameters, this is an issue that needs to be addressed, 
separately from the effect that the inverse of a noisy but 
unbia sed correlation matrix is not unbiased l|Hartlap et al.l 
I2007I ). 



APPENDIX C: SOURCE REDSHIFT 
DISTRIBUTION 

Using lens galaxies at higher redshifts, the mean lens- 
ing efficiency (D\ B /D B ) becomes more sensitive to the 
adopted redshift distribution of the sources, P(z so ). There- 
fore, we have updated the method for determining P(z so ). 
van Uitert et all i|201ll ) used the photometric reds hift cat- 
alogu es of the CFHTLS "Deep Survey" fields (jllbert et all 
2006) and selected all objects in the range 22 < r' < 24 
that satisfied the selection cuts as described in the release 
notes that accompanied the catalogues, i.e. only objects with 
reliable photometry in all the bands, that were observed 
in unmasked regions and with a best-fit template number 
< 54. Since the main interest there was to determine the 
redshift distribution rather than to select galaxies with re- 
liable photometric redshifts, galaxies in the redshift range 
0.05 < Zphot < 2.0 were selected instead of 0.2 < z Dhot < 1.5 
where the redshifts were deemed reliable. Ivan Uitert et all 
(2011) did not account for the scatter of the photometric 



redshifts, nor for the fraction of outliers. Also, they did not 
account for the fact that bright sources have a larger weight 
in the lensing measurements than faint ones. 

To increase the precision of the lensing efficiencies for 
higher lens redshifts, in this paper, we use the photo- 
metric redshift ca talogue from the 2 deg 2 COSMOS field 
|llbert et alj|2009l ) instead. The photometry in 30 bands re- 
sults in photometric redshifts that are both more accurate 
than those from the CFHTLS, and also more reliable up 
to higher redshifts. Using the overlap with the CFHTLS- 
D2 catalogue, kindly provided by H. Hildebrandt, we deter- 
mined the conversion between the r + -band from the COS- 
MOS catalogues, and the r'-band from the CFHTLS. Us- 
ing this conversion we selected source galaxies in COSMOS 
based on their r + magnitudes corresponding to a selection 
of 22 < r' < 24. 

When integrating D\ s /D s over the P(z so ), we have to 
account for the fact that bright galaxies have a larger weight 
in our lensing measurements than faint ones. For this pur- 
pose, we determined the average lensing weight of the source 
galaxies in the RCS2 in narrow r'-band magnitude bins, find- 
ing that on average the sources with r' ~ 22 have a weight 
that is twice that of r' ~ 24 source galaxies. We used the 
conversion between the r + - and r'-band to compute the cor- 
responding weight of each galaxy in the COSMOS catalogue, 
and used that weight to determine the weighted mean lens- 
ing efficiency. 

To account for the outliers, we assigned a new red- 
shift to a random fraction of the galaxies equal to the out- 
lier fraction. The new redshift was drawn from the pho- 
tometric redshift distribution of the sources, and replaced 
the catalogue value when it fulfilled the outlier criterion 
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Table CI. Lensing efficiencies determined from the photometric 
redshift distributions of source galaxies 



^icns {Du/D a ) (D le /D s ) (CFHT) (D ls /D s ) (z phot < 2) 



(1) (2) (3) (4) 



0.1 0.774 0.777 (1.00) 0.768 

0.2 0.602 0.586 (1.03) 0.588 

0.3 0.461 0.440 (1.05) 0.443 

0.4 0.350 0.328 (1.07) 0.329 

0.5 0.264 0.240 (1.10) 0.241 

0.6 0.194 0.172 (1.13) 0.172 

0.7 0.141 0.118 (1.19) 0.118 

0.8 0.103 0.079 (1.30) 0.081 

0.9 0.076 0.053 (1.43) 0.055 

1.0 0.057 0.036 (1.58) 0.037 



(1) lens redshift; (2) lensing efficiency det ermined using the pho- 
tometric redshift catalogues of COSMOS dllbert et aljboOStl ; (3) 
lensing efficiency determined using the photo metric redshift cat - 
alogues of the CFHTLS "Deep Survey" fields jllbert et alj|2006h . 
restricted to source galaxies in the range 2 p hot < 2. The brack- 
eted values show the ratio between column 2 and 3; (4) lensing 
efficiency determined using the photometric redshift catalogues of 
COSMOS, restricting the source galaxies in the range z p hot < 2. 

|.z r andom - «phot|/(l + Random) > 0.15. The outlier fraction 
depends on the brightness; we adopted a value of 0.7% for 
galaxies wit h i + < 23, and 15.3 % for galaxies with i + > 23, 
as quoted in lllbert et all l|2009l ). We created 16 realizations 
of the photometric redshift catalogues, each with a different 
randomly assigned set of outliers, and adopted the mean as 
our new lensing efficiencies. The scatter between the differ- 
ent realizations is small, and can safely be ignored compared 
to the statistical errors of the lensing analysis. We show the 
mean lensing efficiency at 10 lens redshifts in the second 
column of Table [C] 

We ignored the impact of scatter of the photometric red- 
shifts with respect to the spectroscopic redshifts. The effect 
of scatter is that it moves galaxies in redshift from where 
their abundance is large to where it is small. To estimate 
the impact that might have on (D\ s /D s ), we additionally 
scattered each photometric redshift by randomly drawing a 
value from a Gaussian, whose widt h depends on the g alax- 
ies' i-band magnitude, as quoted in lllbert et al.l (|2009l ). We 
multiplied that random value with 1 + z p hot and added it 
to ZphoJEJ We created 16 new realizations, and determined 
the mean lensing efficiency. We found that the impact is less 
than a percent at all lens redshifts, and can therefore be 
safely ignored. 

To see how the lensing efficiencies compare to those 
computed using the CFHTLS "Deep" catalogues, we ap- 
plied the same procedure to compute the average lensing 
efficiencies. We accounted for outliers by adopting the out- 
lier fractions as a function of i'-band magnitude from Ilbert 
et al. (2006), and applied the same weight as a function of 



r'-band magnitude. However, we only selected galaxies with 
Zphot < 2. Again, we created 16 realizations, and determined 
the mean. We show the resulting values of {D\ B /D S ) in the 
third column of Table [Cl At low redshifts, the resulting lens- 
ing efficiencies only differ by a few percent compared to the 
ones based on the COSMOS catalogue. However, we find 
that if the lens redshift increases, the {Di a /D s } from COS- 
MOS becomes increasingly larger. To demonstrate that this 
difference is due to source galaxies at z p hot > 2, we repeated 
the calculation using the COSMOS photometric redshift cat- 
alogue, but now restricting the analysis to z p hot < 2. We 
show the resulting lensing efficiencies in the fourth column 
of Table [C] We find that the lensing efficiencies agree very 
well with those based on the CFHTLS "Deep" catalogues. 
The difference is at most 4% over the entire redshift range 
that we probed. 

In previous work where we us e d the photometric red- 
shift catalogues from lllbert et al.l (|2006l ) to compute the 
lensing efficiencies, we focused at galaxies at low redshifts. 
Hence the lensing efficiencies that we used there were of suf- 
ficient accuracy. However, for galaxies at redshifts z > 0.5, 
our results show that is it important to include source galax- 
ies at Zphot > 2 in the computation of {Di„/D a ). 

Note that we have ignored cosmic variance. However, 
we find very similar lensing efficiencies using the COSMOS 
and CFHTLS "Deep" photometric redshift catalogues when 
we restrict the galaxies to z p hot < 2. This suggests that 
cosmic variance does not have a large impact on the lensing 
efficiencies that we use. 



11 Formally, lllbert et ail ||2009| 1 quote the scatter on Az/(1 + 
Zspec) with Az = Zphot — Zapoc, so we should have multiplied the 
random value with 1 + z spe c rather than 1 + z p i lot . However, we 
expect the difference to be minor. 
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